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Abstract 

We study the effect that conjugation-mediated Horizontal Gene Transfer (HGT) has on the mutation-selection balance of a 
population in a static environment. We consider a model whereby a population of unicellular organisms, capable of 
conjugation, comes to mutation-selection balance in the presence of an antibiotic, which induces a first-order death rate 
constant k d for genomes that are not resistant. We explicitly take into consideration the repression/de-repression dynamics 
of the conjugative plasmid, and assume that a de-repressed plasmid remains temporarily de-repressed after copying itself 
into another cell. We assume that both repression and de-repression are characterized by first-order rate constants £ + _and 
k- + , respectively. We find that conjugation has a deleterious effect on the mean fitness of the population, suggesting that 
HGT does not provide a selective advantage in a static environment, but is rather only useful for adapting to new 
environments. This effect can be ameliorated by repression, suggesting that while HGT is not necessarily advantageous for a 
population in a static environment, its deleterious effect on the mean fitness can be negated via repression. Therefore, it is 
likely that HGT is much more advantageous in a dynamic landscape. Furthermore, in the limiting case of a vanishing 
spontaneous de-repression rate constant, we find that the fraction of conjugators in the population undergoes a phase 
transition as a function of population density. Below a critical population density, the fraction of conjugators is zero, while 
above this critical population density the fraction of conjugators rises continuously to one. Our model for conjugation- 
mediated HGT is related to models of infectious disease dynamics, where the conjugators play the role of the infected (I) 
class, and the non-conjugators play the role of the susceptible (S) class. 
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Introduction 

Horizontal Gene Transfer (HGT) is defined as the transfer or 
exchange of genetic information between two organisms, that does 
not involve the transmission of genetic information from parent to 
daughter as a result of replication, a process known as vertical 
transmission of genetic information [1,2]. Since it was initially 
characterized in the 1940s [3], it has gradually become evident 
that HGT is one of the major sources of genomic change in 
prokaryotes [2]. As a result of its recognized importance for 
shaping the genomes of prokaryotes, a complete theory describing 
the evolutionary dynamics of prokaryotes must take HGT into 
account [2]. Further, the fact that HGT plays a major role in 
shaping prokaryotic genomes means that HGT must be consid- 
ered when reconstructing phylogenetic trees [1], which of course 
complicates any effort to understand the evolutionary relationships 
amongst bacterial strains, and to trace the evolutionary histories of 
various genes. The existence of HGT also has consequences for 
public health, since HGT is believed to be primarily responsible 
for the rapid spread of antibiotic drug resistance in bacterial 
populations [2,4-8]. 



Despite the large influence that HGT has on the evolutionary 
dynamics of prokaryotic populations [2], with a few notable 
exceptions [9-12], there has been relatively little theoretical work 
attempting to incorporate HGT into evolutionary models. Indeed, 
it is safe to say that the bulk of the work in evolutionary biology has 
focused on point-mutations as the primary source of genetic 
variation [13]. There are two main reasons for this: First of all, the 
simplest evolutionary models are the ones that consider point- 
mutations as the only source of genetic variation. Yet even these 
models can exhibit fairly non-trivial behavior, and there is a wide 
range of additional effects that may be considered within the 
point-mutation framework (e.g., conservative versus semiconser- 
vative replication, fitness landscapes, polysomic genomes, genetic 
instability, genetic repair and mutators, finite-size effects [9,13- 
21]). Thus, already at the simplest level of point-mutations, there is 
a fairly large class of evolutionary models that may be analyzed. 
Second, for quite some time in biology it was believed that point- 
mutations are indeed the primary source of genetic change in 
organisms. This assumption underlies the concept of the so-called 
"tree of life" that was used to organize phylogenetic data. With the 
realization of the importance of HGT, there have been suggestions 
to replace the tree analogy with that of a web or a ring [1,22—25]. 
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Bacterial genomes typically consist of a single, large, circular 
chromosome, in addition to several smaller, circular chromosomes 
known as plasmids. Bacterial plasmids can move between bacteria, 
often between bacteria of different strains, via a process known as 
conjugation. [26-29] . If these plasmids contain genes for an adaptive 
trait such as antibiotic drug resistance, then conjugation can 
significantly speed up the rate at which such an adaptive trait 
spreads throughout a population. [30-33]. Conjugation is believed 
to be one of the most important forms of HGT, so that both 
phenomenological and quantitative models of processes such as 
the emergence of antibiotic drug resistance will need to properly 
account for this behavior. 

The best characterized bacterial conjugation system is the FV 
F system [28]. Here, a bacterium containing what is termed an 
F-plasmid fuses with a bacterium lacking the F-plasmid. The 
bacterium containing the F-plasmid is termed an F + bacterium, 
while the bacterium that does not contain this plasmid is termed 
an F bacterium. When an bacterium meets an F bacterium, 
it transfers one of the strands of the F-plasmid to the F bacterium 
via a pilus. Once a strand of the s F-plasmid has been transferred 
from the F + bacterium to the F bacterium, a copy of the plasmid 
in both cells is produced by daughter strand synthesis using the 
DNA template strands. The F bacterium then becomes an F + 
bacterium that transcribes its own pilus, and is able to transfer the 
F plasmid to other bacteria in the population (Figure 1). 

The F + /F system is in some ways atypical for bacterial 
conjugation systems, due to its conjugative machinery being 
chronically in a de-repressed state, while in nature most 
conjugative plasmid systems are repressed via the expression of 
T4SS modules [34,35]. Nevertheless, the F-plasmid system is still 
one of the best characterized bacterial conjugation systems, and it 
can be considered as broadly representative of all known bacterial 
conjugation systems. Therefore, it makes sense to base mathe- 
matical models of conjugation-mediated HGT on the F-plasmid 
system. 

In recent work [36], the authors developed mathematical 
models describing the role that conjugation has on the mean 
fitness of unicellular populations in static environments, and 
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Figure 1. The conjugation process for F+ bacterium with F- 
bacterium as detailed in our model: 1. A de-repressed bacterium 
comes within proximity of an F- bacterium/repressed F+ bacterium. 2. 
The F+ bacterium's pili interacts with the F- bacterium's membrane. 3. 
The F+ plasmid releases a strand of its DNA to the F- bacterium via 
conjugation machinery. 4. A new de-repressed F+ bacterium is formed. 
doi:1 0.1 371 /journal, pone.0096839.g001 



found that conjugation-mediated HGT does not confer a 
selective advantage to the bacterial population (indeed conjuga- 
tion-mediated HGT was found to have a slightly deleterious 
effect on the mean fitness of a population). However, in this 
work, it was assumed that the conjugation machinery is 
permanently de-repressed. Here, we extend our previous work 
by introducing the repression/derepression mechanism into our 
original model [13] of a bacterial population that transfers 
genetic information among its members both vertically and 
horizontally (Figure 2). 

Results 

1. Definition of the Model 

Our model assumes a unicellular population of asexually 
reproducing bacteria, where the bacterial genome consists of 
two chromosomes. The first chromosome is the main circular 
bacterial chromosome that transfers its genetic information via 
vertical transmission (i.e. replication from one generation to the 
next). The second chromosome is taken to be the much smaller 
conjugative plasmid, which is assumed to consist of two regions: 
(1) A "conjugation" region, comprising the various genes 
necessary for bacterial conjugation. (2) An "antibiotic resis- 
tance" region, comprising the genes necessary for conferring 
resistance to a given antibiotic. Since each chromosome is 
assumed to be a double-stranded, semiconservatively replicat- 
ing DNA molecule, we define probabilities for the error-free 
replication of a template strand from each of the various 
regions of the bacterial genome. To this end, we let p v denote 
the probability that a template strand from the main bacterial 
chromosome produces a daughter chromosome without intro- 
ducing any new point-mutations. We define p c and p r 
analogously for the portions of the bacterial plasmid coding 
for conjugation and antibiotic drug resistance. We assume that 
there are master sequences for the main bacterial genome, the 
conjugation and antibiotic drug resistance regions of the 
plasmid. We assume that any mutation to these master 
sequences renders the sequences non-functional. We further 
assume that the bacterial genome controls viability, and so a 
bacterium with a master, or wild-type, copy of its main 
chromosome, has a normalized first-order growth rate constant 
of 1, while a bacterium with a mutated main chromosome has a 
first-order growth rate constant of 0. The assumption that a 
single mutation renders the master sequence non-functional is 
known as the single-fitness-peak approximation [13]. Although 
it is an oversimplification of the true dependence of function- 
ality on the gene sequence, it nevertheless captures the fact that 
only a small fraction of sequences are likely to yield a gene that 
encodes the desired function. As a result, this approximation 
has been known to give quantitative results in certain cases 
[37]. 

Given its simplicity, it therefore makes sense to start with this 
assumption before moving on to more complicated fitness models. 
We assume that a viable bacterium whose resistance region of the 
plasmid corresponds to the wild-type is unaffected by the presence 
of an antibiotic, while a viable bacterium with a mutated portion 
of this resistance region is killed by the antibiotic at a rate 
characterized by a first-order rate constant Kj). 

Bacterial conjugation is modeled as a second-order rate process 
characterized by a rate constant y. We also assume that, as the 
size of the bacterial population changes, the system volume 
changes so as to maintain a constant population density p while 
assuming that p remains constant. It should be noted that in this 
context, the only bacteria whose population we are tracking are 
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Figure 2. Examples for vertical replication and horizontal gene transfer within the framework of our model and the consequential 
propogation of mutations throughout the population: 1 . Mutation-free vertical replication of the de-repressed, fully viable subpopulation that 
is notated as + . 2. Vertical replication followed by a mutation to the conjugation and antibiotic resistance genes in the plasmid. 3. Mutation-free 
bacterial conjugation between an F+ and an F- bacterium that leads to the formation of two antibiotic resistant, conjugation capable, derepressed 
bacteria. 4. Horizontal transfer between an F+ and F- bacterium followed by a mutation to the antibiotic resistance gene of both conjugants, leading 
to no drug resistance for these bacteria. /? w respectively notate the probabilities for the bacterium replicating its viability, conjugation, and 
resistance portions of the genome in an error-free manner while yp is the second order reaction parameter of the conjugation rate. 
doi:10.1371/journal.pone.0096839.g002 



the viable ones. The idea behind this assumption is that the 
bacteria each take up a certain amount of space, and so, as the 
bacterial population grows, the total system volume must increase 
as well. Bacterial conjugation can only occur between a viable 
bacterium whose conjugation genes are de-repressed, and a 
viable bacterium whose conjugation genes are repressed. The 
repression/de-repression dynamics is assumed to be character- 
ized by first-order rate constants £ + _and k- + , respectively, 
corresponding to transitions from the de-repressed to the 
repressed state, and from the repressed state to the de-repressed 
state, respectively. However, for the purposes of this paper, we 
will assume that k-+ =0, corresponding to the limit of weak de- 
repression, which is believed to characterize the repression/de- 
repression dynamics of actual systems [34,38]. 

During bacterial conjugation, we assume that a template strand 
of the plasmid from the de-repressed bacterium, which is labelled 
the "+" bacterium, is transferred to the repressed, or " — ", 
bacterium, and that daughter strand synthesis occurs in both 
bacteria to reconstitute the plasmids. We assume that the plasmid 
transferred from the "+" bacterium replaces the plasmid in the " — 



" bacterium, and inherits the de-repressed state of the parent 
plasmid. The assumption of plasmid replacement is a simplifica- 
tion that will need to be re-examined in future research, where we 
anticipate developing more accurate models that allow for variable 
plasmid numbers in the bacterial cell. The basis for this 
assumption derives from the observation that plasmids of similar 
compatibility classes cannot co-exist in the same cell [39,40], and 
that bacteria can control the number of plasmids in the cell [40— 
42]. While these observations do not directly translate into a 
plasmid replacement model, we believe that they do provide some 
justification for this assumption. 

Finally, we assume that replication errors are due to mismatches 
during daughter strand synthesis that are sub-sequentially fixed in 
the genome. We let s denote the mismatch probability per base- 
pair, so that 8/ 2 is the probability of making a mismatch, which is 
fixed as a mutation in the genome. If L v ,L c ,L r denote the lengths 
of the viability, conjugation, and resistance portions of the 
genome, respectively, then we have that p v = (1 — e/2) Lv , 
p c = (\-&/2) Lc <mdp r = (l-s/2) L \ 
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Now, following the standard procedure in population genetics 
and quasispecies models [13,19], we may allow the genome length 

L = L V + L c + L r to become very large while holding \i = sL ^ + = [2p v p c p r — l—k + -+ yp(2p c p r — \)x c _ — fc(t)]x\ + 

constant. Defining a v = L v /L,a c = L c /L,a r = L r /L, we then ob- 

cix 

tal11 ' = [2PvPcPr ~ 1 -ypx\ ~ K(t)]x+ + + k+ _ X+ + 

Pv,c,r = e-«^ 2 (1) ^±=, =[2p vPc -l- KD -k + .+ YP (2p c -l) X c _ -K(t)}<- + 

2Pc(l-Pr)xX+(P» + yPX°- ) 

2. Evolutionary Dynamics Equations d ^ = [2 PvPc -l - KD -ypx c + -ic(0]x;_ + 

Here, we define x An to be the fraction of the population "t 

consisting of viable bacteria such that the conjugator region is of k + -X+_ + 2p v p c (\ —p r )x+ + 

type A, and the resistance region is of type □, where a type of "+" ^ 

means that the given region is of the wild-type, while a type of " — " — = [2p v p r — 1 — ypx c + — k(t)]x- + + 
means that the given region differs from the wild-type. When 

A = +, corresponding to a conjugator, then °= ±, corresponding 2p v (l —p c )p r (x+ + 

to the genes for conjugation being de-repressed (° = +) or repressed ^ x i + ^ c 

(° = — ). Since we are only tracking viable bacteria, we should note ^ ^ c P r i x + + X - 

that xl, + x7 , +xt_ +x~_ =x dx-- r . _ ... 

•wit j £• c 11 ir • r • i i i_' • i j = P^v-l-K^-ypX 0 , -K(/)X— + 

We define x + to be the total traction ol viable bacteria that are at 

conjugators, so that x c + =x+ + +x| _ , and we define x c _ to be the 2 p v (\ -p c )(\ -p r )(x% + +x~ + ) 
total fraction of viable bacteria that are non-conjugators, so that + _ 

x c _=x~ + +x~_+x- + +x--. Therefore, allowing us to sum- +2p v (\ -p c )(x + _ + x + _) + 2p v (\ -p r )x_ + + 

marize all of the processes that are relevant for each subpopulation yp[2(l —p c )]x c [(1 —p r )x+ + +x^_ } 
in our model into a generalized list of balance equations which are 
described by Eq. (2): 



(3) 



error-free semi-conservative replication (vertical gene transfer) : 

+ 2 PvPjPk-^ + + 
X ±± } X ±± +X ±± 

non-plasmid-restricted mutation (as a result of vertical gene transfer) : 

+ 2 P v PjPk + ^ + ^ 
X ±± ^ X ±± +X ±± 

error-free bacterial conjugation (horizontal gene trasnfer) : 

YpVPjPk ~ 1) 

x\ + x c _ > x\ + x c + (2) 

plasmid-restricted mutation (as a result of bacterial conjugation) : 

x c + +x c _ > x c _ +x c _ 

repression : 




death by antibiotics : 



x J _ — > dead bacterium 



We note that pj denotes two possible values for the conjugation 
gene - p c when the gene is copied or conjugated in an error-free 
manner, or 1 — p c when the gene is copied incorrectly. In a similar 
manner, p^ denotes two possible values for the antibiotics 
resistance gene - p r or 1 —p r . 

Putting everything together, we obtain that the evolutionary 
dynamics of the bacterial population is governed by the following 
system of equations: 



By summarizing all of these equations we get that 
K(t) = 2p v — 1 — tcd(x+_ +x~_ is the mean fitness of 

the population, or the first-order growth rate constant of the 
population as a whole. Finally, we'd also like to note that the form 
of the equations detailed by Eq. (3) is identical to what would be 
obtained if we made a more realistic assumption that the 
population was growing in a chemostat. [13]. 

3. Mutation-selection Balance 

We solve our system of equations for the mutation-selection 
balance for the case when ji is small enough so to allow us to avoid 
the error-catastrophe region. This allows us to apply a first order 
taylor expansion to Eq. (1), which gives us simpler terms for p v ,p c 
and p r which are p v = 1 —a v /n/2;p c = 1 — a c p/2 and p r = 1 —a r jn/2 
(the probabilities for an error-free replication of bacterium's 
chromosome, conjugation gene and resistance gene respectively). 
Our first goal is to determine the initial conditions which will allow 
our system to converge while maintaining physical context. We 
proceed by considering a preliminary condition for the antibiotics- 
related death/Co which must be greater then the value of zero 

while denoting Xq + + ,Xq + + ,Xq + _ ,Xq + _ ,Xo, - + ,Xo, as the 

fi = 0 values of x+ + ,x~ + ,x+ _,x~ _,X- + ,x , respectively. We 

show in sub-section A from the Materials and Methods that this 
condition brings forth that there are no non-resistant organisms in 

the system when ji = 0 so that Xq + _ =Xq + _ =Xo, =0. From 

this point, we separate our analysis into two regions based on the 

differences between the values ofypand The first region 

obeys the condition of yp <k + _ and we have proved that for this 
region Xq + + = Xq + + = Oand Xo, - + =1 (sub-section B from the 
Materials and methods). The second region obeys the condition of 
yp>k + - and we proved that for this region 
Xq + + = 1 —k + _ / (yp),XQ + + =k+-/ (yp),Xo, - + = 0 (sub-section 
G of the Materials and methods). We applied these to Eq. (3) and 
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found two separate solutions for and consequent- 

ly were able to compute the terms for the population mean fitness 
for each of these regions (sub-section D in Materials and methods). 
For yp<k + _ , we have Xq + + =Xq + + =0,and Xq,~ + = 1 while 
x have the following values: 



= 0 



x~_ = 0 



2(1 - Pr ) 
K D 



(4) 



Therefore, in selection-mutation balance their population mean 
fitness is: 



K = 1 — (0t v + 0C r )/H 



(5) 



The second solution is for yp>k + - when 
Xq ++ = 1 — k + - /yp,XQ ++ =k + -/yp and Xo,-+=0 when we 
obtain, 

<--^o-^x. + *-> 

c _ _ 2{\-p r ) k+-{yp-k+-){\+k+-) + k+-K D (6) 
' L 7P(i<d+7P-I< ) 



K D 



x__ =0 



from which we derive the selection-mutation balance mean 
fitness of: 



, = ,-i..„,,.-^>.^-H '-^ | - 1 ^ +t ;; H - 



(7) 



In order to analyze the effect of the conjugation rate yp on k we 
took the first derivative with respect to yp, and observed that k is a 
decreasing function of yp, going from 1 — (oc v + oc r )fi for yp = k + _ 
to 1 — /n[ot v + ot r (\ +k+ _)] as yp ->oo (Figure 3). 

We proceeded with a similar analysis for the repression rate 
k + - and calculated the first and second derivatives of the 
population mean fitness/c with respect to k + -. We found that 
k + - =yp + KD — \/ Kj)(yp-\-Kj)) is a local minimum of k 
(Figure 4). We can then take the limit of »0, and find that 

k=\ — /i(oc v + a r ) and that k is a decreasing function for 
k + - <k + -^ C rit. We treat the limit of »yp of k similarly 

and find that at the range of k + -, C rit <k + _ <yp,k + _ reaches the 
exact same value of k as in the case of k + - ->0 (Figure 4). 

Discussion and Conclusions 

We have solved our model for the limiting case of weak de- 
repression k-+ = 0, and found a phase transition between a non- 
conjugator and a conjugator regime (Figures 3, 5, and 6). In our 
previous work [36], we found a similar phase transition in a 
system, where all of its population was de-repressed. Such a 
system, can also be considered as the strong de-repression limit for 



our current model (k+ _ ->oo, k + _ =0). Nevertheless, a phase 
transition that is related to the conjugation rate in both cases was 
derived independently for each one of them. Therefore, we 
conjecture that such a phase transition is ubiquitous whenever the 
rate of repression is smaller than the rate of conjugation (i.e. 
0<£ + _<yp). 

At first glance, the phase transition that we found for our 
current model may seem somewhat counter-intuitive, since we 
would expect no conjugators for all values of yp due to the effect of 
repression. Nevertheless, we can explain the appearance of 
conjugators from the transient de-repression of the " — " cell upon 
conjugation. When we examine the trends of the population 
fractions given in the model with respect to this transition, we 
observe that when yp is sufficiently large compared to the 
characteristic rate of repression, the transient de-repression results 
in a chain reaction of conjugation-induced de-repressions 
throughout the population, leading to a stable, and positive, 
population fraction of conjugators. In other words, once the 
conjugation rate yp crosses a certain threshold, the transient 
production of conjugators becomes sufficiently large enough to 
make the steady-state with no conjugators an unstable one. Any 
perturbation away from this steady-state takes the system to a 
steady-state where conjugators are present in the population. 
Interestingly, we noted from Eq. (14) that the threshold value for 
yp (the rate of conjugation) increases with k + - (the rate of 
repression) and decreases with Kd (the antibiotics-induced death 
rate), sub-section E in the Materials and Methods). This is 
consistent with the idea, that in order to establish a stable 
conjugators population, the rate of conjugation has to be slower 
than the rate of death caused by the antibiotics, as to allow the 
removal of the non conjugators population via mutation (due to 
signal fitness peak landscape, this would effectively lead to the 
removal of the non conjugators) while at the same time, the rate of 
conjugation has to be significant enough to avoid its effective 
cessation by repression. 

When we introduced bacterial conjugation into our population, 
the mean fitness of the population in our model has decreased. 
This is in general agreement with previously published experi- 
mental results [43], which attribute this fitness cost to its 
maintenance by the bacteria. A similar effect was observed in 
our previous work on the subject [36], although in that specific 
case the mean fitness has decreased to a minimum for a critical 
value of the conjugation rate yp, and asymptotically increased to 
the value of k as yp—» go (with respect to its value for k(yp crit )). 

However, in our current model we have found that, the mean 
fitness decreases monotonically from 1— (a v + a r )/i to 
1 — (a v + (X r )p — k + -(X r p without any observed increase of the 
mean fitness at the limit of yp— ► oo. We believe that this may imply 
on an HGT- driven phase transition that resembles the error- 
catastrophe phenomenon in viruses that could lead to a critical 
decrease of the population's mean fitness [44], However, the scope 
of such a phenomenon needs to be further investigated in further 
theoretical as well as experimental studies. 

Furthermore, we noticed that the difference in the values of 

/cincreases with the value of the repression rate k^ This 

dependency, implies on the complicated relationship among the 
parameters yp, k + _ and /cwhich we have further investigated by 
differentiating k with respect to k + _ as well as by yp. From these 
derivations we have found that as the conjugation rate is increased, 
the value of the mean fitness decreases. However, this decrease can 
be effectively nullified if the rate of repression is high enough. This 
is the result of a local minimum for k + _ (the repression rate) that 
occurs at a critical value, which we obtained from differentiating 
Eq. (12) with respect to k+- (the derivation of this value can be 
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Figure 3. A plot of k (the population mean fitness) versus yp (the second order reaction parameter of the conjugation rate). 

Parameter values area v = 0.8,a c = a r = 0.1,K;.D = 1.0,fc + _ = 10.0,// =0.1 {oc c /oc r /oc v are the ratio of the genome which is dedicated to conjugation/ 
resistance/organism viability respectively, kd is the first-order death rate constant of the population that is caused by antibiotics, k + - is the first 
order rate constant which describes the conjugators transition from the de-repressed to the repressed state and describes the probability of 
mutation for every replication of the entire genome). The solid line represents the analytical expression for k obtained in the small-ja limit, while the 
pluses represent the values for k obtained numerically. The numerical values for k were obtained using fourth-order Runge-Kutta integration of the 
evolutionary dynamics equations. For yp<k + -, we used the initial condition x_+ = 1.0,x+ + = x~ + , = x+_ = x~_ =x__ =0.0, while for yp>k + - 
we used the initial condition x+ + = 1.0— k+-/(yp),x~ + = k + - /{yp),x\_ = x~_ =x_+ =x__ =0.0. We iterated for 100,000 time steps of size 
0.001. 

doi:10.1371/journal.pone.0096839.g003 




Figure 4. A graphical solution for the population mean fitness K versus k (the first order rate constant corresponding to the 
transition from the de-repressed to the repressed state of the conjugators). The values of k were calculated from Eq. (12) with the 
parameter values: a v = 0.8, a c = a r = 0.1, kd = 1.0, yp= 10.0, ^ = 0.1. 
doi:10.1371/journal.pone.0096839.g004 
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100 



Figure 5. A plot of population fraction of conjugators x c + + versus the conjugation reaction parameter yp. Parameter values 

are the same as those for figure 3, including the parameters for the Runge-Kutta integration. The solid line represents the zeroth-order analytical 
expression for x c + , while the pluses represent the values obtained numerically. 
doi:1 0.1 371 /journal. pone.0096839.g005 



found in sub-section F in the Materials and Methods and is 
illustrated by the graphical solution of Teas a function of k+ _ in 
Figure 4). For k+- <k + -^ cr i t the value of yp is big enough to 
lower the population mean fitness while for k + _ >k + -^ cr i t the 
rate of transition from the de-repressed state to the repressed state 
is sufficiently large to decrease the difference in the mean fitness 
caused directly by conjugation. We should note that both the 



k + -^>yp and »0 limits of k give the same value of 

K= 1 — /z(a v + a r ). Therefore, once the repression rate is greater 
than the value ofk + -^ cr i t = yp + ?Cd — yj Kj)(yp-\-Kj)), the decrease 
to the mean fitness caused by conjugation is mitigated to the point 
that it is effectively cancelled by repression. 

Our study illustrates the critical role that repression is likely to 
have on regulating the process of conjugation and how it can 



08 - 



0.6 - 




0.4 - 



0.2 - 



100 



YP 



Figure 6. A plot of the population fraction of the non conjugating, antibiotic resistant bacteria sub-population t versus the 
conjugation reaction parameter. The parameter values are the same as those for figure 3, including the parameters for the Runge-Kutta 
integration. 

doi:10.1371/journal.pone.0096839.g006 
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negate the deleterious effect on the fitness of a population. 
Therefore, while conjugation is disadvantageous in a static 
environment, its negative effect on the fitness of a population is 
likely to become negligible if that population can sufficiently 
repress the plasmid and consequently, restore itself to its original 
fitness. This would suggest that the role of the bacteria in the 
actual host-parasite dynamics between the bacteria and the 
plasmid are much more complicated. While it is known that 
bacteria can obtain various advantageous traits from plasmids 
(including antibiotic resistance [3], degradation of unusual 
substances such as xylene, toluene, napthtalene and phenanthrene 
[45], resistance to heavy metals [46,47], and translation of Colicins 
[48]), they are also required to pay a fitness cost for hosting the 
plasmid in a static landscape. Previously, it was suggested that 
repression could reduce the energy costs for the host bacteria [49] . 
However, experimental results suggested that the amelioration of 
the fitness cost due to plasmid hosting was possible through 
evolution of a bacterial population [50] . Based on our results, we 
propose an alternative mechanism that allows bacteria to 
ameliorate this fitness cost by repressing their propensity towards 
conjugation and not only by vertical evolution. Interestingly, this 
strategy is likely to allow bacterial populations a significantly faster 
restoration pathway for their fitness loss (due to the maintenance 
costs of antibiotic drug resistance) then vertical gene transfer (i.e. 
evolution). 

Our current model bears similarities to the Susceptible-Infected- 
Susceptible (SIS) [51-53] class of infectious disease models, where 
the " — " cells play the role of the Susceptibles (S) and the "+" cells 
play the role of the Infected (I). If the contact frequency between 
the S and the I population classes is sufficiently low, so that the Fs 
revert back to the S state before having an opportunity to infect at 
least one additional S, then the fraction of I's in the population 
drops to a negligible fraction of the population. Above a critical 
contact frequency, however, the fraction of I's rises to a finite 
positive fraction of the population. Our model is in broad 
agreement with the behavior of actual conjugative plasmid 
systems, where it was observed that conjugative plasmids can 
spread in an infectious manner above a critical population density 
[38,54,55]. Typically this decrease in the population mean fitness 
is explained by the presence of Male Specific Phages (MSPs) [56] , 
but in our model this factor is not accounted for, resulting in a 
decrease of the mean fitness that is a direct result of mutations that 
the plasmids go through when being replicated within the host 
bacteria. 

We would like to also note recent studies that have suggested 
bacterial conjugation as a platform for the development of 
conjugation-based antibiotics [57]. In light of these studies, our 
results may have implications about the pharmaceutical limits that 
need to be considered for such antibiotic agents. 

Finally, we believe that developing mathematical models that 
can correlate known data is important, as it indicates that the basic 
mathematical framework is reasonable, and may therefore serve as 
a foundation for building more sophisticated quantitative and 
predictive models. For future work, we will need to take into 
account more biological features of conjugation, such as regulation 
of plasmid copy number, MSPs and plasmid compatibility classes. 
Furthermore, we will need to move away from static landscapes 
and consider the role that conjugation has in allowing a population 
to adapt to new environments. 



Materials and Methods 

A. Proof That There Are No Non-resistant Organisms for 

H = 0,kd >0 

If we define x r + to be the total fraction of resistant organisms, 
and x r _ to be the total fraction of non-resistant organisms, then we 

have x r + =x| + and x r _ =x|_ +x~_ +x We 

then have that fi = 0, 



dt 

dx~ + 

dt 
dx\ _ 

dt 
dx~ _ 

dt 
dx- + 

dt 
dx-- 

dt 



■■[l—k+- + ypx c _ —K(t)]xX + 
--[\-ypx\ -K(t)]x- + +k+-xl + 

■ [1 — kd — k + _ + ypx c _ — K{t)]x~X _ 

= [1 -K D -ypx\ —K(t)]x+_+k+-xX- 
--[\-ypx\ -K(t)]x- + 

■ [1 —Kj)—ypX C + — K(t)]X-- 



(8) 



If we assume that the system converges to a stable steady-state, 
then we must have that K~(t)>\ — k + - +ypx c _ ,1 — ypx\ . If, for 
Kj)>0 there are non-resistant organisms in the population, then 

either x + _ ,x~ _ on are strictly positive. If x+ _ >0, then we 

have k= 1 — Kd — k + _ +ypx c _ < 1 — k + - +ypx c _ =><=. There- 
fore, x^_=0, so either x~^_ or x are >0, which gives 

K= 1 — Kd — ypx c _ < 1 — ypx c _ =^><^.Thus, there are no non-resis- 
tant organisms in the population at steady-state for ji = 0. 



B. Proof That xi 



v o,+ - 



= 0,x 0 ,- + = 1 for yp<k^ 



Let us suppose that x+ + > 0 for someyp <k + _ Then the first 
line in Eq. (3) gives k = 2p v p c p r — 1 — k + _ + yp(2p c p r — \)x c _ . Now, 
assuming that our system converges to a stable steady-state, we also 
have from the the fifth line in Eq. (3) made > 2p v p r — 1 — ypx\ , and 
so yp-k + - >yp[x\ +(2p c p r -\)x c _]-k + ->2p v p r (\-p c )> 
0^><^. Therefore, x^ + =0 for yp<k + -,/n, which implies that 



x, 



o,+ + 



= 0. 



Suppose then that, x~ + >0for yp<k + - ,{i>0. Then from Eq. 
(3) we have/c = 2p v p c p r — 1 — ypx c + < 2p v p r — 1 — ypx c + =><=, 
hence x~ + = 0 for yp <k + _ ,fi>0. 

Taking the limit as we obtain x$ ++ =0for yp<k + -. 

Now, when fi = 0, we have Xq + + +Xq + + +Xo,- + = 1, and so for 
we obtain Xq _ + = 1 . 



C. Proof That x^ 



= \ -k A 



o,++ - A -^+-/(yp)v*o,++ = k +-/(yp\ 

xo,-+ =0 for yp>k + - 

Given some yp>k + _ , suppose that Xq + + =0. Then since that 
are no non-resistant organisms for/i = 0, we have that Xq + = 
x C q ++ =0 and Xq _ = 1 . This implies thatx > 1 — k + _ + yp> 
1^><^, and so we must have that x^ ++ >0. This gives, 1 = K = 
l-k+- +yp(^oT++ +*o,- + )=>*£++ +*o,-+ = k+-/(yp), and 

so x o,++ = l ~ k +-/(yp). 
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Now, the fifth line of Eq. (3) gives, 

0= -ypx 0 + + + x 0 ,- + 



(9) 



from which it follows that Xo,-+=0, and so 

+ =k+-/(yp). 
D. Derivation of the Small-|i Expressions for 

For small JI, we may write the first-order expressions for 
dx^_/dt and dx--/dt. We obtain, 



0= [-K D -k+- +yp(x^ ++ +x 0 ,_ + )]x+_ 
+ 2(1 -p r )x£ + + (l + yp(x 0 ~; + + +* 0 ,- + )) 
0= -[K D + ypx^ + + ]x~_ + 2(l-/? r )x 0 ~ + + 

0= -[/cz> + ypx 0 + + + ]x__ + 2(l-/? r )x 0 ,- + 



(10) 



For yp<k + -,we have x^ + + = Xq + + = 0,and Xo,-+ = l, 
which implies = x~_ =0, and. 

X-- =x r _ =2(1 —p r )/Kj)i thus K = 2p v — 1 —Kj)X r _ =2p v — 1 — 
2(l-p r ) = l-(a v + a r ) i u. 

For yp>k + -\ve have x 0 h ++ = 1 — k + -/yp,XQ + + = k+-/yp, 
and Xq 5 -+ =0. This gives, 



2(1-^) 


k + - 


yp 


k +- 

2(1 -i?r) 


a-'- 


1 



x__ =0, 



Equating this term to zero allows us to find (yp) CJ 



dK K D k + - (2yp + K D -k + -) 
= — pa r — — — — = 0 



diyp) 

-+2yp + K D -k + - = 0 

^(yp)crit=\( k +-- K D) 



(yp) 2 (yp + K D -k+-) 2 



(14) 



F. Derivation of k + - yCrit 

In order to find the critical value of ^ + _we take the derivative 
of Eq. (12) with respect to £ + _6>and get the following expression: 



dk, 

dk 



Ik d f , 
— = — jl-Kav + a, 



yp i _ _|_ _/p 



a,-p(yp + k 0 ) [(yp) 2 + yp(K D -2k+-)-k+- (2k d -£+_)] 
yp(yp + K D -k+-) 2 



(15) 



By equating this expression to zero we get the following 
quadratic equation for k + _ : 



oc r p(yp + fcz>) [(yp) 2 + yp(K D - 2k + _)-k+- (2k d -k+-)]= 0 



(yp) +yp(K D -2k + -)-k + -(2K D -k + -) = Q 
ypK D - 2ypk + _ - 2^/;+ _+/:+_ 2 + (yp) 2 = 0 
k + - 2 -k + - (2yp + 2tcd) + yp(yp + k d ) = 0, 



(16) 



which allows us to find the values for k+- that solve the 
aforementioned equation: 



and so the mean fitness at steady-state is given by, 



k=\-h\ a v + a r [(l-^— )(l+fc + _)+^e /V_ vo " 



(12) 



2yp + 2k d ± yj(-2yp- 2k d ) 2 -4-1 -ypiyp + k^) 
2l : 



]) 2yp + 2k d ± y/4yp 2 + Syp-^ + 4k d 2 -4yp 2 - Ayp-K D _ ( 17 ) 



yp^K D ±y^(yp~^KD) 



E. Derivation of (yp) crit for yp>k + _ 

In order to find the critical value of yp, we compute the 
derivative of Eq. (12) with respect toypoand derive the following 
expression: 



Once we take into account the singularity of K at 1 — + ^ 
we are left only with the following: 



k+ - tCri t = yp + K D - \/K D (yp + K D ) 



(18) 



dK d 
d(yp) diyp) 



1 — fi((X v + oc r 



K D k + - (2yp + K D -k + -) 
(ypf{yp + K D -k + -f 



(13) 
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Note that one can compute that. 

d 2 k 2(x r 'KD'fi(yp + K£>) 2 
dk\_ yp(yp + K D -k + -) 3 

d 2 k[yp H~ Kp \J K D (yp + k d )] _ 2u r -K D -ii(yp + K D ) 2 
dk + _ 2 yp(yp + K D -k + _) 3 

_ 20L r 'K D 'li(yp + K D ) 2 _ 

ypiyp + k d - [yp + k d - \jK D (yp + K D )] ) 3 

2a r 'K D -[i(yp + K D ) 2 2a r -K D -fi(yp + K D ) 2 

ypiyp + K D -yp-K D + ^K D (yp + K D )f yp(^K D (yp + K D )f 
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If we take into account the positivity of yp,Kj),tt r ,n, we see that 
2.,K D .,(yp + K D y _ ^ so ^ cn . cal 

ypW K D(yp+KD)) 



yp + Kp — \J Kd(yP H~ Kd)is indeed a (local) minimum value of K. 
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